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Abstract — This paper addresses a simplified frequency offset estimator for multiple -input-multiple^^lput 
(MIMO) orthogonal frequency division multiplexing (OFDM) systems over frequency-seleo ffi^ e fading 
channels. By exploiting the good correlation property of the training sequences, which are coja^riiifed from 
the Chu sequence, carrier frequency offset (CFO) estimation is obtained through factor d^o&position for 
the derivative of the cost function with great complexity reduction. The mean squaiWerroV (MSE) of the 
CFO estimation is derived to optimize the key parameter of the simplified estimatoA^Well as to evaluate 
the estimator performance. Simulation results confirm the good performance of JfeeVyJuning-assisted CFO 
estimator. XV^? 

Index Terms — Frequency offset estimation, frequency-selective fading idafl^ieis, low complexity, multiple- 
input- multiple-output (MIMO) orthogonal frequency-division mul^^^^OFDM). 

I. Introduction /^V-/ 

Orthogonal frequency division multiplexing (OFDM) transmission is receiving increasing attention in 
recent years due to its robustness to frequency-selective foi^wj^and its sub carrier-wise adaptability. On the 
other hand, multiple-input multiple-output (MIMO) sVfcAris attract considerable interest due to the higher 
capacity and spectral efficiency that they can providau^omparison with single-input single-output (SISO) 
systems. Accordingly, MIMO-OFDM has emerget^^ 1 * strong candidate for beyond third generation (B3G) 
mobile wide-band communications [i].It isja^ll known that SISO-OFDM is highly sensitive to carrier 
frequency offset (CFO), and accurate estim^ra^jand compensation of CFO is very important [2]. A number 
of approaches have dealt with CFO estii|f9&^& in a SISO-OFDM setup [3], [4], [2], [5], [6], [7]. According to 
whether the CFO estimators use tr#n^g*§cquences or not, they can be classified as blind ones [3] [4] and 
training-based ones [2], [5], [6],J7]?S^rfclar to SISO-OFDM, MIMO-OFDM is also very sensitive to CFO. 
Moreover, for MIMO-OFDM, t|p£te exists multi-antenna interference (MAI) between the received signals 
from different transmit ant yy ta^The MAI makes CFO estimation more difficult, and a careful training 
sequence design is requiratf^jjtraining-based CFO estimation. However, unlike SISO-OFDM, only a few 
works on CFO estimatjofl^j^ MIMO-OFDM have appeared in the literature. In [8], a blind kurtosis-based 
CFO estimator for m/j€>-OFDM was developed. For training-based CFO estimators, the overviews 
concerning the «efcsj^ry changes to the training sequences and the corresponding CFO estimators when 
extending SISQ^l^lvI to MIMO-OFDM were provided in [9], [10]. However, with the provided training 
sequences in^L satisfactory CFO estimation performance cannot be achieved. With the training sequences 
in [10], tJjjaS^iining period grows linearly with the number of transmit antennas, which results in an 
increa^?^overhead. In [11], a white sequence based maximum likelihood (ML) CFO estimator was 
'for MIMO, while a hopping pilot based CFO estimator was proposed for MIMO-OFDM in [12]. 
taal calculations of the CFO estimators in [11] [12] require a large point discrete Fourier transform 
(DFT) operation and a time consuming line search over a large set of frequency grids, which make the 
estimation computationally prohibitive. To reduce complexity, computationally efficient CFO estimation 
was introduced in [13] by exploiting proper approximations. However, the CFO estimator in [13] is only 
applied to flat-fading MIMO channels. 

When training sequence design for CFO estimation is concerned, it has received relatively little attention. It 
was investigated for single antenna systems in [14], where a white sequence was found to minimize the 
worst-case asymptotic Cramer-Rao bound (CRB). Recently, an improved training sequence and structure 
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design was developed in [15] by exploiting the CRB and received training signal statistics. In [16], training 
sequences were designed for CFO estimation in MIMO systems using a channel-independent CRB. In [17], 
the effect of CFO was incorporated into the mean square error (MSE) optimal training sequence designs for 
MIMO-OFDM channel estimation in [18]. Note that optimal training sequence design for MIMO-OFDM 
CFO estimation in frequency selective fading channels is still an open problem. 

In this paper, by further investigating the above search-free approaches, a simplified CFO estimator is 
developed for multiple input- multiple - output (MIMO) OFDM systems over frequency selective fading 
channels. With the aid of the training sequences generated from the Chu sequence [9], we prc»po£^o 
estimate the CFO via a simple polynomial factor. Thus, the complicated polynomial rooting operf^i^ is 
avoided. Correspondingly, the CFO estimator can be implemented via simple additions and multoglic^nons. 
To optimize the key parameter of the simplified CFO estimator as well as to evaluate Jby^ltimator 
performance, the mean square error (MSE) of the CFO estimation is derived ^S>^ 

Notations: Upper (lower) bold-face letters are used for matrices (column vectors). S^^stripts , T and H 
denote conjugate, transpose and Hermitian transpose, respectively. (-)P denotes tJjeVg^due of the number 
within the brackets, denote the floor, Euclidean norm-square, expectation^^d Kronecker product 
operators, respectively, sign(-) denotes the signum function and sign(o) = 1 is *^Vied. [x]m denotes the m- 
th entry of a column vector x. x(m) denotes the mcyclic-down-shift veraQM^f^:. diag/x/ denotes a diagonal 
matrix with the elements of x on its diagonal. [X]m,n denote the (m^j^Jntry of a matrix X. FN and IN 
denote the N x N unitary DFT matrix and the N x N identity mato|^^pectively. E k N denotes the k-th 
column vector of IN. 1Q (oQ) and oPx Q denote the Qx 1 all-one Cafl^z^o) vector and P x Q all-zero matrix, 
respectively. JQ denotes the Q x Q exchange matrix with oneswi its anti-diagonal and zeros elsewhere. 
Unless otherwise stated, o < |i < Nt - 1 and o < v < Nr - 1 are^jumed 

II. Signall 

Consider a MIMO OFDM system with Nt transn%jraJ«ennas, Nr receive antennas, and N sub carriers. The 
training sequences for CFO estimation are the^ameas in [6] and [7]. Let s denote a length-P Chu sequence 
[9]. Then, the P x 1 pilot sequence vector atJffl^ph transmit antenna is generated from s as follows: 

K = y/Q/NtF P a^ M where M ^^efine 




^9 — l e JV> e JV ' • • • • e N ^> 




Then, the N x 1 training s|fJ^Sice vector at the (ith transmit antenna is constructed as follows: 



Let y v denote the i^^^Jceived vector at the vth receive antenna after cyclic prefix (CP) removal. Let h(v,|_i) 
denote the L x i%^mj^l impulse response 




a Q - L and i p = y iff /i = p. For 
"er to {t^Q 1 as the Chu-sequence-based train- 



vector, with L being the maximum channel length. Assume that L is shorter than the length of CP Ng. Let ~e 
denote the frequency offset normalized by the subcarrier frequency spacing. Define 
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T T T T "\T 

vi>y[ vi y,v r -ij 



(hMy,(hwy t 
...,(fc<^)) T ,...(/ l ("' JV '- , )) T ] T 

^O-^l h„ ,...,h Nr _i\ 

D„(i) =diag {[l.e* 2 ^/* eiS^N-D/N]^! 

Then, the cascaded received vector y over the Nr receive antennas can be written as [6], [7] 

y= JNe*' eN »' N {I Nr ®[D N {f)S}}h + w (l) iSO^ 

Where ^> % 

■•^e^- 1 ®er Nt _ 1 |{l Nf s{F J v[/,,0 J , x(N _ L) ] T |} >S>r 
and w is an NrN x 1 vector of uncorrelated complex Gaussian noij^ymples with a mean of zero and an 
equal variance of 02W ^ti 

HI. CFO Estimator for Mtt^TOFDM Systems 

By exploiting the periodicity property of CBTS, y can^fc^acked into the Q x NrP matrix 

Y = [To, Yi, . . . , Y v, . . . ,YNr- 1] with its element grvfen by 
X*]i* = [({e~ Nr ) T 0 /iv)y],P + p 

define # 



B(i) = ffc 0 , 61 • • f^ft* ■ - il • 

Then, Yean bej^^i^ed in the following equivalent form [6], [7]: 
Y = B(iY!jQW (2) 

where X 1 ^^ 

>/. [»M>, . . . . . . 

x Fjdiag{^}O^F N [lL,0 Lx(N . L) ] T h^ 

and W is the Q x NrP matrix generated from w in the same way as Y According to the multivariate 
statistical theory, the log-likelihood function of Y conditioned on B(e) and X, with e denoting a candidate 
CFO, can be obtained as follows: 
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\np(Y\B(e),X) = -<7-*TV{[y - B(e)X] [V - B(e)X] M } . 



(3) 



Exploit the condition i\x = i\x_ iff [_i = Then, after some straightforward manipulations, we can obtain the 
reformulated log-likelihood function conditioned on e as follows: 



\np(Y\e) = Tr[B H (e)R YY B(e)] 



(4) 



h Direct grid searching from 
where R YY =YY (4) 



Yields the ML estimate; however, this approach is computationally quite expensive. To efficien^Jbmpute 
the CFO, we will subsequently propose a simplified CFO estimator for MIMO OFDM sysfrefesjD efine z = 
^ne/Q^ z ^ _ gjzmii/Q^ an( j = ^ z Z q f z q_ Then, by exploiting the F^mrl^in p 

the log-likelihood function in (4) R Y Yean be transformed into the following equiva"" 



Where 




mman property of, 
rtn: 



qth element given by [rjV =q[* Ryy ]hj • It ( 
corresponds to the s^i^Mhiion of the qth upper di 
'especf 

f'{z) = z- 1 |c T j ^ G 6(z) 0 q| 



Where c is a Q x 1 vector with its qth 
definition that the qth element of c 

YY Taking the first-order derivative of f[z) with respeclfc^: yields 



can be seen from its 
diagonal elements of R 



where q = [o, 1, . . . , q, . . . , Q - ffiVBy letting the derivative of the log-likelihood function/1 (z) be zero, the 
solutions for all local mimraffl^rnaxima can be obtained. Put these solutions back into the original log- 
likelihood function /(z), |9lrtwect the maximum by comparing all the solutions obtained in the previous 
stage. The improved hjm*S&FO estimator that exploits the above mathematical rule has been addressed for 
single-antenna OFQM^^tems in [3]. Although the search-free approach has a relatively lower complexity, 
it still requires aNcAnj^icated polynomial rooting operation, which is hard to implement in practical OFDM 
systems. With/r^^id of the CBTS training sequences, we will show in the following that the polynomial 
rooting opejjInSn can be avoided for training-aided CFO estimation in MIMO OFDM systems. Assume that 
P > L, tftiXrannel taps remain constant during the training period, and the channel energy is mainly 
conx^rj^^ea in the first M taps, with M <L. Then 




with 1 < l < Q - 1, and the parameter 1 denotes the index of the upper diagonal of . From (7), it R Y Y 
immediately follows that/ (z) can be decomposed as follows: 
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(8) 



It follows from (8) and (9) that z = ~ z is one of the roots of both/ (~ z) = o and zQ - k(l) = o. Unlike f_{~ z) 
= o, the roots of zQ - k(l) = o can be calculated without the polynomial rooting operation. Therefore, by 
solving the simple polynomial equation zQ - k(l) = o, the CFO estimate can be efficiently obtained as 
follows: 

A- 

m^ilxity of the 

whose main 
simplified CFO 



£= argmax {f(z)\z = e^ € ^} 



(10) 



Where eq = arg/K( 1)^/(271) + q - Q/ 2. It can be calculated that the main computat 
simplified CFO estimator is 4NrNQ+ 8Q2. Compared with the CFO estimator in [ 
computational complexity is 4NrN log2 N + 9Q3 + 64/^Q - 1)3, the complexity 



estimator is generally lower. Furthermore, since the polynomial rooting operatio 
CFO estimator can be implemented via simple additions and multiplication 
practical OFDM systems. Note that 1 is a key parameter for the proposed CF' 
the following how to determine the optimal l. 





the simplified 
is more suitable for 
imator. We will show in 



IV. Simulation Results^ 

Numerical results are provided to verify the analytical results as^^l as to evaluate the performance of the 
proposed CFO estimator. The considered MIMO OFDM s^\rn^has a bandwidth of 20 MHz and a carrier 
frequency of 5 GHz with N = 1024 and Ng = 80. Each of trA^mnels is with six independent Rayleigh fading 
taps, whose relative average powers and propagation dfcAs are {o y - 0.9,- 4.9,- 8.0,- 7.8,- 23. 9 } dB and /o, 
4, 16, 24, 46, 74^ samples, respectively. The other paxamAers are given as follows: P = 64, Q = 16, Nt = 3, Nr = 
2, and ~ e G (-Q/ 2,Q/2). Figs. 1 and 2 present the^jSrrf the proposed CFO estimator as a function of l with 
{i\i}Nt- 1 \i=o = {3, 5, 11} and fe, 7, 14}, resp^ively. The solid and dotted curves are the results from the 
analysis and Monte Carlo simulations, resf^o^vely. It can be observed that the results from the analysis 
agree quite well with those from the siraffllttSsns, except when the actual MSE of the estimate is very large. 
It can also be observed that MSEf" #iitrtoe/es its minimum for 1 = 6, 8, 10 with {i\i}Nt- 1 (i=o = fe, 5, 11} and 
for l = 7, 9 with {i\i}Nt- 1 (i=o = 7, Tlfc^Vrhese observations imply that we can obtain the optimum value of 
the parameter 1 from the analyt^l^esults after {i\i}Nt- 1 \i=o is determined 




Fig. 1. MSE of the proposed CFO estimator as a function of 1 with 
iNt i={3, 5, 11}. 
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Fig. 2. MSE of the proposed CFO estimator as a function of l 
with? Ntl ={3,7,14}' 
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Fig. l CFO estimation performance for diffqrcH^m 

We resort to Monte Carlo simulatjtrt/iSwts evaluation. It can be observed that the performance of the 
proposed estimator with CBTS k far^r^er than that in [10] and slightly worse than that in [6] and [7], and 
its performance also approachi^the EMCB, which verifies its high estimation accuracy. It can also be 
observed that the performanj«3£*ne proposed CFO estimator with CBTS is far better than that with RS, 
which should be attribute^Aftl^ good correlation property of CBTS. 

V. Conclusion 

In this paper,^j^fi9ve presented a low-complexity CFO estimator for MIMO OFDM systems with the 
training seqi<Ace5 generated from the Chu sequence. TheMSE of the CFO estimation has been developed 
to evalua^^e estimator performance and to optimize the key parameter. By exploiting the optimized 
paramff£kn%m the estimation MSE, our CFO estimator with CBTS yields good performance 

Or 
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